Stretching force dependent transitions in single stranded DNA 
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Mechanical properties of DNA, in particular their stretch dependent extension and their loop formation char- 
acteristics, have been recognized as an effective probe for understanding the possible biochemical role played by 
them in a living cell. Single stranded DNA (ssDNA), which, till recently was presumed to be an simple flexible 
polymer continues to spring surprises. Synthetic ssDNA, like polydA (polydeoxyadenosines) has revealed an 
intriguing force-extension (FX) behavior exhibiting two plateaus, absent in polydT (polydeoxythymidines) for 
example. Loop closing time in polydA had also been found to scale exponentially with inverse temperature, 
unexpected from generic models of homopolymers. Here we present a new model for polydA which incoipo- 
rates both a helix-coil transition and a over-stretching transition, accounting for the two plateaus. Using transfer 
matrix calculation and Monte-Carlo simulation we show that the model reproduces different sets of experi- 
mental observations, quantitatively. It also predicts interesting reentrant behavior in the temperature-extension 
characteristics of polydA, which is yet to be verified experimentally. 

PACS numbers: 87.10.Pq,87.15.La,05.70.Fh 
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In order to understand how the mechanical properties of 
DNA and RNA influence biological processes like transcrip- 
tion and translation in living cells and viruses, these biopoly- 
mers are stretched in vitro to study their nonlinear elasticity 
and their internal structure. ssDNA, despite receiving less at- 
tention compared to double stranded DNA (dsDNA) ^ ^, 
has recently attracted a lot of interest. Smith et al. pl] showed 
that the FX diagram for wild type A-phase ssDNA can be de- 
scribed by FJC model only at low force (< IQpN). To ex- 
plain the behavior at higher force they used a modified freely 
jointed chain (mPJC) model yj, [SO with stretch dependent 
kuhn length. Subsequently, it was discovered that synthetic 
ssDNA have interesting sequence dependent properties. For 
example, polydA was found to have a higher bending rigid- 
ity than polydT and as a consequence polydA takes relatively 
longer time to form a loop [J4|]. Further, the loop forma- 
tion time varies exponentially with inverse temperature which 
can neither be explained by FJC nor worm like chain (WLC) 
model, which are generic models for flexible and semiflexible 
polymers, respectively. This behavior was attributed to stro_n 
stacking interaction among pyrimidine bases in polydA 
|6[]. Subsequently, FX characteristics of polydA revealed that 
it undergoes two successive transitions under external stretch- 
ing ItZi Isll generating two distinct plateaus in the force exten- 
sion (FX) curve. The first one at ^ 23pN force was proposed 
to be a helix to coil transition in which the inter-base stack- 
ings are broken and helical polydA segments transform to a 
polydA coil. The second transition at ^ llipN force was at- 
tributed to the over-stretching of the constituent bases. It was 
conjectured 131 17|] that over-stretching results from the con- 
formational change of the sugar molecules from C3' -endo to 
C2' -endo pucker conformation. 

Zimm-Bragg model jgl] which was originally proposed to 
explain temperature driven helix-coil transition in proteins, 
has been used to study the force driven helix-coil transition 
in polydA. Theoretical models have used mean field approxi- 
mation inn and exact evaluation of partition function [12y to 



explain the first plateau experimentally seen jlOll in the FX 
diagram of polydA at low forces (< 60pN). But the sec- 
ond plateau, involving the over-stretching transition is beyond 
the scope of these models. Overstretched in dsDNA has been 
studied using Ising like two state models by various groups 
i 131 1 1411 . Here we propose a model which quantitatively re- 
produce both the force driven behavior as well as zero force 
conformational fluctuations like the loop formation time, as 
observed experimentall y. D ouble plateau behavior has also 
been addressed before lllSi Il6|] . theoretically, albeit using a 



lattice model. But lattice models generically underestimate 
entropic effects and also quantitative comparisons to experi- 
ments were not possible. 

Model: We model polydA as a chain of connected seg- 
ments, each representing a nucleotide. Length of a segment 
(bond) represent the phosphate to phosphate distance in the 
ssDNA backbone. The bond length can be Ih, Ic or Is depend- 
ing on whether the segment is in the heUx {Ih = Q.yinm 
llOD), coil {Ic = 0.59nm ||2[]), or overstretched-coil state 
(Is = 0.7nm Im)- While the value of Ic has been most 
well documented ||2|] in the literature, Ih is determined lliol] 
by noticing the 1.6 times increase in contour length of polydA 
upon helix to coil transition. Here Ih of course is the pro- 
jected length of the helical contour on the central axis of the 
helix. Further, Is = O.lnm has been inferred from the maxi- 
mum extension reached by the polydA chain under very high 
force (- eOOpiV) in Ref|0, i]. Incidentafly, Is = 0.7nm 
also matches with the distance between the consecutive phos- 
phates when the deoxyribofuranose ring is in the C2' -endo 
pucker conformation. In our model, state of a segment is 
characterized by (yu, S), where yu can take values or 1, cor- 
responding to coil or helix state, respectively. Once in the 
coil state ^ = 0, there are two possible states: the normal 
coil or the overstretched coil, corresponding to bonds lengths 
Ic and Is- These two states are represented by S* = —1 and 
S = +1, respectively. Hence, there are only three possible 
states (1, —1), (0, —1) and (0, 1), corresponding to helix, coil 



and overstretched-coil. The Hamiltonian of this three state 
model of polydA (omitting the external force) is 

N 

Ho =^(2(l-^,)/i,+iAw + /x,+iA/ 

~JSA+i{l - mO(1 - M.+i) + hS,{l - Ai.)) (1) 

This Hamiltonian incorporates two transitions, actually cross- 
overs, since the model is one dimensional. First part of the 
Hamiltonian involving Aw and A/, is the simpler version 
of the original Zimm-Bragg model H, forwarded by the au- 
thors themselves and later used by Tamashiro et al. 1I12I1 in the 
context of forced DNA. In this model the necessary hydrogen 
bonding, that is required for the formation of helical domain, 
takes place between adjacent bases (segments), instead of the 
i — th and (z + 4) — th bases as in the original Zimm-Bragg 
model ||9[]. The weight of the configurations cc, he, ch, hh, 
in this model, are given by 1,1, as and s respectively. Here 
a = e-2A«;/fcBT ^jjjj ^ ^ g-AZ/feflT^ ^j^gj-g 2Aw is the inter- 
facial energy between the helical and the coil domains and A/ 
is the difference of energy between the helix and the coil state. 
The asymmetry between ch and he arises because a segment 
(its oxygen atom) can engage either its right neighbor (its hy- 
drogen atom) or its left neighbor to make a hydrogen bond 
(for details see Ref|l9|]). Second part of the Hamiltonian is an 
Ising Hamiltonian, which will be invoked, when the segments 
are in the coiled or overstretched coil state i.e., when ji — Q. J 
measures the correlation energy between an adjacent coil and 
overstretched coil states. 2h is the energy difference between 
a coil and a overstretched coil state. The parameters a and J 
are often called the cooperativity parameters of the respective 
transitions. 

We considered the ssDNA chain to be semi-flexible, in 
which the helical domains have very large persistence length, 
where as the coiled and overstretched-coil domains have small 
persistence length. Discretized worm-like-chain hamiltonian 
for the system is 

JV-l 
PHbend = ^ l—fiifii+1 + (1 - Hi){l - l-k+l) X 



(1 - 5.)(1 - 5,+i) + y(1 + S,)il + 5,+i)]} 



(1 -cos9i^.i+i). 



(2) 



Here, 0i,i+i is the angle between the bond vectors ti and 
<i+i, where ti = Ri+i — Ri and R.; is position vector of 
i — th monomer. The hamiltonian has been constructed in a 
way such that different bending rigidities are associated when 
neighboring bonds are of identical type, i.e., hh, cc or ss. It 
amounts to assuming that helical, coil and stretched-coil do- 
mains of the polymer behave like worm-like-chain, ah , dc 
and as are the respective persistence lengths. We choose 
ah ~ 12nm, a^ = 1.5nm and Ug = l.bnm, i.e., a relatively 
large persistence length for the helix (still much smaller than 
that of dsDNA, about bOnm). Theoretical models (l^ [H] 



assume ah to be infinite. Our assumption that even ssDNA 
coils have a small persistence length is in agreement with 
Seol et al. 's JITIl FX data onpolyU where they obtained 
Uc ^ Inm. Also Smith et al |3y reported a^ ^ 1.5nm for 
wild type A-phase ssDNA, using mFJC model. 

For calculation purpose, we substitute inter-bond angles in 
terms of bond vectors: ^(1 — cos9i^i+i) = f (ii+i — h)^. 
Finally, including external stretch F, the total hamiltonian is 



Hi 



total 



— Ho + Hhend — F-(RAr — Rq) 



(3) 



The force dependent term above can be expressed 
as F. J2i=i ^i^i ™ terms of the bond vectors and 
their respective lengths li, where li could be lh,lc 
or Is depending on the internal state (/i, S) of the 
bond. Finally, the partition function is Z{F) = 
/ ^tidii E{^.,s.} exp hl^iHo + H^end - FR,)], where 
Rx = (Rw — Ro)-a; is projection of the end-to-end distance 
along the force F = Fx. Using transfer matrix technique we 
can write Z{F) for the forced chain as. 
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/^i,5i,ii;/j.iv,5jv,ijv 



{^ll,Sl,tl\r'-'\^lN,SN,tN) (4) 



Here the transfer matrix T is an 3mn x 3mn matrix where 
the internal state space is 3 dimensional, corresponding to 
the states (1, — 1), (0, — 1), (0, 1) and the orientation space 
ti{(j), 9) has been discretized into m. x n bins. In EqU Z{F) 
is obtained as a weighted sum over all the matrix elements of 
y w-i We cannot exploit the simplifications normally arising 
from periodic boundary condition because, here, the polymer 
has one of its ends fixed and the other end stretched by a force. 
The details of this calculation can be found in Ref 1 1 811 . After 
calculating Z{F) we can compute (R^) as a function of force 
using (Rx) = 4 ^i*^^ . This is shown in Fig. 1. 

Force-extension behavior : In general the stretching force 
tends to align the chain along the force (x), at the cost of en- 
tropy. It would also favor the individual bonds to have their 
highest possible bond lengths i.e., O.lnm in order to maxi- 
mize the FRx term. This is achieved only at very high force 
when the other terms in the Hamiltonian give in to the force 
term. But at low and intermediate forces the other terms com- 
pete. Although helical segments are favored over coil seg- 
ments (due to s), at low force and at room temperature, en- 
tropy has substantial contribution. As a result all the seg- 
ments are not aligned along the force. That is why the ex- 
tension per base, at very low force is about 0.2jim and not 
0.37nm (see Fig.l). But bending of the helical domain is 
disfavored by its large persistence length. Therefore, at low 
force, the entropy induced bends help some helical segments 
to convert to coil, which have higher internal energy but lower 
bending energy compared to helical segments. As the force 
rises, the FRx term enforces alignment as well as stretching 
of the bonds, leading to stacking-unstacking transition which 
adds extra length (that was previously curled up in helix) to 
Rx- The abrupt nature of the helix-coil transition results from 
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FIG. 1: (color online) Force-Extension curve obtained from our 
model (solid line), using transfer matrix method, compared with the 
experimental result of Chen et al. [SJI (circle). Such a double plateau 
feature was first reported by Ke et al and their data is very close 
to that of ReflSQ shown here. The inset shows how fraction of the 
different species, namely, helix(circle), coil(square) and stretched 
coil (triangle) portions change with increasing force, obtained from 
Monte-Carlo simulation of our model. The intersection points be- 
tween the helix and coil fractions mark the position of the first plateau 
(~ 23pN) and that between the coil and overstretched-coil fractions 
mark the second plateau ('^ 114pA'') in the FX-plot. Parameters used 
for our calculation are 2Aw = GMbT, A/ = -4.93fci3r, J = 
OAiksT, h = l.bksT, where T is the room temperature. We used 
two different chain lengths, Ai" = 33 and 65, to check A'^ indepen- 
dence of this plot. 



coupling term /^ti/ii+i in the Hamiltonian. The width of the 
plateau, i.e., how much length is released at the transition, de- 
pends on how many segments are converted from the helical 
to the coil state. This comes from the abruptness of the jump, 
shown in the inset. We compute the area under the FX curve: 
/ F.dx ^ GksT. Given that 2h = 'iksT is the free-energy 
change for the second transition, the helix-coil transition ac- 
counts for a free-energy change of about SfegT. 

Given that polydT and polyU do not have any significant 
stacking interaction, they might be expected to follow the 
A/ = Aw — limit of our model. But polydT,polyU 
also have qualitatively different overstretching behavior than 
polydA; they overstretch very slowly with force (unlike the 
transition in polydA). PolyU has been already explained by a 
modified WLC model lIlTIl and we checked that even polydT 
can also be explained by the same model but with different 
parameters. 

Loop formation kinetics : Now we discuss the zero force 
conformational fluctuations resulting from our model. We 
focus on the observations made by Goddard et al 141] on 
loop formation properties of short (N = 8 — 30 bases) 
polydA and polydT chains. They had attached complemen- 
tary base sequences TTGCC and AACGG at the two ends 
of a poly d A/poly dT strand and attached flurophore (F) and 
quencher (Q) molecules at the ends. This design aimed to de- 
tect the formation of hairpin loops by zipping of complemen- 



tary base-pairs at the ends. Such a process resulted in quench- 
ing of fluorescent intensity of F by Q. They found that, for 
a given chain length, a polydA chain took longer time than a 
polydT chain to form a loop. This indicates greater bending 
stiffness for polydA chains, resulting from stacking interac- 
tion between Adenine bases. More intriguing was the result 
that, loop closing time In(Tc) oc (3 for polydA and a nearly flat 
temperature dependence for polydT. 

This process can be approximately described by a two state 
system: a chain in a open state (o) or a chain with a closed 
loop(c). At equilibrium, the interconversion o ^ c obeys 
a detail balance condition: pckc^o ~ Poko-^c |6|], where 
Pc/o are the equilibrium densities and kc^o, ko^c are the 
conversion rates. Loop closing and opening times, Tc and 
To, are inverses of the respective rates. Assuming a small 
interaction radius a between the chain ends, approximately 
Pc = ^a^PN{R = 0), where PNiR) is the probability 
of finding the chain ends at a separation R. But since pos- 
sible number of open configurations far outweighs number 
of closed configurations i.e., po » pc, we can approximate 
Po = 1 — Pc — 1 la] . Further, Tq is determined by the high 
energy barrier of the five bases long zipping strand, which is 
independent of the chain length N, and hence is a constant. 
Thus we arrive at t^ oc p~^, which we can compute from our 
model, as a function of chain length (N) and temperature, at 
zero force. In Fig. 2 we plot Tc/N^^^ versus inverse tempera- 
ture, obtained from our model and compare it with the exper- 
imental data of Goddard et al ]|4[]. The rationale for rescaling 
Tc by N^^'^ is to partially nullify the strong A^ dependence in 
Tc- Although Pc oc N~'^/'^ only for an FJC model (at large 
N) and not for an WLC model, nevertheless it turns out to be 
useful in approximately collapsing both our simulation data 
and Goddard et al's experimental data in a narrow range of N , 
near the room temperature. 

Temperature dependence of extension : Change of exten- 
sion with temperature has been shown to give interesting be- 
havior for wild type ssDNA 111911 . Fig. 3 shows the analo- 
gous property resulting from our model of polydA, at fixed 
force, computed using transfer matrix method. In the absence 
of any transitions polymer extension is expected to decrease 
with temperature because of entropic elasticity, as it happens 
in rubber, for example. In case of wild type ssDNA hairpin 
loops can form which modifies the extension-temperature be- 
havior in non-trivial ways. For polydA although loops can- 
not form in the absence of complementary bases, existence 
of two transitions (helix-coil and overstretching) makes it be- 
have in an interesting way (see Fig3 A), to the extent that weak, 
nonmonotic, re-entrant behavior can be observed (Fig3B,C). 
These can be understood qualitatively. Two important points 
to remember here are, a) helix-coil transition can also be af- 
fected by raising temperature, and b) entropic fluctuations 
are enhanced at high temperatures which smoothens the force 
driven transitions (data not shown). Fig3A shows that at very 
low force like 5pN, indeed extension weakly decreases with 
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FIG. 2: (color online) Semi-log plots of the scaled loop closing 
time Tc (of polydA) versus inverse temperature, computed for var- 
ious chain lengths A'^ = 10 (circle), 20 (diamond) and 30 (triangle), 
by Monte-Carlo simulation of our model. In a narrow window of 
temperature lOOO/T £ [3.1,3.45], the scaling of Tc has been ex- 
perimentally shown (inset-A) to be approximately, Arrhenius type 
(i.e., Tc (X exp(e/fc_BT)) by Goddard et al 0]. Their data for chain 
lengths A'^ = 12 (circle), 21 (diamond) and 30 (triangle), are plotted 
after rescaling with A'^'' "). Both the Arrhenius aspect as well as the 
quantitative value of the exponent e are closely reproduced by our 
model. The dashed lines in both the main plot and the inset have the 
same exponent e = 13. 4fcs To, where To = 2981-4^ is the room tem- 
perature. Inset-B schematically shows transition between open and 
closed chain conformations. 



FIG. 3: (color online) Extension(E) per base (in nni) versus temper- 
ature, at a fixed force, plotted in (A) for different values of forces: 
5, 10, 15, 20, 25, 30, 60, 100 and 150piV, from bottom to top. These 
are obtained from our model using transfer matrix method. Effect 
of the first transition (near 25pN) is reflected in the relatively large 
jump in the extension. Near the overstretching transition the jump 
is not so pronounced due to the relatively small plateau width in the 
corresponding FX diagram. The variation of extension with tem- 
perature changes qualitatively with force. The extension weakly de- 
creases with temperature at low force, increases with temperature at 
intermediate force (< 30pA'^), shows interesting reentrant behavior 
at higher force (30 — lOOpA'^) and finally again decreases with tem- 
perature at very high force. The reentrant behavior i.e., extension 
initially increasing with temperature and later at higher temperature, 
decreasing with temperature, is zoomed in (B) and (C). The minute 
change in the FX curves with temperature is shown in (D). Dashed 
lines are guide for the eyes. 



temperature but when force rises ( lOpN onwards) as we move 
close to the helix-coil transition, it is easier to affect the tran- 
sition by raising temperature and helix-coil transition leads 
to rise in the extension. But at forces just beyond the hehx- 
coil transition (30pN onwards but much below lOOpN) the ss- 
DNA cannot access the overstretched bond lengths solely by 
means of thermal fluctuations and looses out to entropic elas- 
ticity showing decrease in extension with temperature. But as 
force approaches lOOpA^ due to its vicinity to the overstretch- 
ing transition the extension again increases first but eventually 
loose out to entropic elasticity at higher temperature, giving 
rise to reentrant behavior Beyond the overstretching transi- 
tion extension again decreases weakly with temperature due 
to entropic elasticity, weakly because at such high force not 
much entropy is left in the almost straight configuration. 

In conclusion, we have proposed a new model for polydA, 
that incorporates two transitions and quantitatively reproduces 
both force-extension characteristics and loop closing statistics 
of such homopolymers. Our model also predicts interesting 
reentrant behavior in the temperature-extension diagram of 
polydA which can be verified experimentally. 

We thank Dibyendu Das for useful comments. 
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